Phase effects in neutrino conversions during a supernova shock wave 
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Neutrinos escaping from a core collapse supernova a few seconds after bounce pass through the 
shock wave, where they may encounter one or more resonances corresponding to Am^tm- The neu- 
trino mass eigenstates in matter may stay coherent between these multiple resonances, giving rise to 
oscillations in the survival probabilities of neutrino species. We provide an analytical approximation 
to these inevitable phase effects, that relates the density profile of the shock wave to the oscillation 
pattern. The phase effects are present only if the multiple resonances encountered by neutrinos 
are semi-adiabatic, which typically happens for 10~^ < sin^ 613 < 10"'^. The observability of these 
oscillations is severely limited by the inability of the detectors to reconstruct the neutrino energy 
faithfully. For typical shock wave profiles, the detection of these phase effects seems rather unlikely. 
However, if the effects are indeed identified in the Ve spectra, they would establish inverted hierarchy 
and a nonzero value of 613. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

The neutrino fluxes emitted from a core collapse super- 
nova (SN) contain information about the primary fluxes 
produced inside the star, the neutrino mixing pattern 
as well as the matter densities that the neutrinos have 
passed through. The high statistics neutrino signal that 
one expects from a future galactic SN needs to be de- 
coded in order to extract this information. 

The neutrinos and antineutrinos produced inside the 
SN pass through the core, mantle and envelope of the 
progenitor star before escaping. They encounter matter 
densities ranging from nuclear densities deep inside the 
core all the way to vanishingly small densities in the in- 
terstellar space. Neutrino masses and mixing angles, and 
hence the extent of their flavor conversions, depend cru- 
cially on the density of surrounding matter, hence it is 
important to study these matter effects in detail. 

The matter effects [l| give rise to MSW resonances [1| 
when the matter density corresponds to 

PR = ±Am'^cos29 MN/{2V2GFYeE) . (1) 

Here Yg is the electron fraction, Mn the average nu- 
cleon mass, and the plus and minus signs correspond 
to neutrinos and antineutrinos respectively. For neutri- 
nos of energy E, resonances are possible at two den- 
sities, the H resonance corresponding to (Ato^,6') w 
(Am^jjjj, 6'i3) and the L resonance corresponding to 
« (Ami,, 012). The energies of SN neutrinos 
are typically in the range 5-50 MeV. For these energies, 
the H resonance takes place around pn ^ 500-5000 g/cc. 
It occurs in neutrinos for normal hierarchy and in an- 
tineutrinos for inverted hierarchy. The L resonance that 
takes place around pL ~ 20-200 g/cc always occurs in 
neutrinos. However, since 9i2 is large, significant flavor 
conversions of antineutrinos also take place at the L res- 
onance. 

The adiabaticities at the H and L resonances deter- 
mine the net neutrino conversion probabilities. The 



L resonance is always adiabatic, given the values of 
AttIq , 012 and the typical density profile of the progenitor 
star around pL- The adiabaticity at the H resonance is 
very sensitive to the value of ^13 and the density profile 
of the star in the resonance region. Indeed the neutrino 
conversion rates are crucially dependent on the value of 
013, and whether the H resonance is in neutrinos or an- 
tineutrinos. The SN signal can therefore be an extremely 
sensitive probe of 0i3 and whether the mass hierarchy is 
normal or inverted [3|. 

In addition to divulging the neutrino mixing scenario, 
SN neutrino fiuxes can also allow us to have a peek at 
the propagation of the shock wave while it is still inside 
the mantle of the star. The violent density fiuctuations 
caused by the SN shock wave can change the adiabatic- 
ity at the H resonance in a time dependent manner, 
thus leaving their imprint on the time dependent neu- 
trino spectra [1, [1,0,13] ■ In particular, the observations 
of the time dependent neutrino spectra can confirm the 
presence of forward as well as reverse shock wave through 
the "double dip" feature [8|, and in addition can track 
the positions of the shocks as they pass through the H 
resonance region. The feasibility of such a trackin g a t a 
water Cherenkov detector has been explored in [ol. llOj. 

Our understanding of the SN explosion mechanism is 
still unsatisfactory [ill [l^ , which makes it very impor- 
tant to extract as much information about the shock wave 
as possible. In this paper, we demonstrate how the neu- 
trinos that pass through the shock wave near the H res- 
onance carry information about the density profile of the 
shock wave. In addition to making the H resonance tem- 
porarily nonadiabatic, the shock wave also forces the neu- 
trinos to encounter multiple H resonances. The relative 
phases that the neutrino mass eigenstates gain between 
two or more of such H resonances manifest themselves as 
oscillations in the neutrino fiavor conversion probabilities 
as a function of the neutrino energy. These oscillations 
are related to the shock wave density profile, and in prin- 
ciple carry information about it. We therefore study the 
relation between the oscillations and the shock wave pro- 
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file in detail. 

Most of the analyses of SN neutrino conversions till 
now calculate the probabilities of conversions between 
neutrino mass eigenstates in matter in each resonance 
region, and combine the results at different resonances, 
assuming them to be independent. This includes the im- 
plicit assumption that the information on the relative 
phase between the neutrino mass eigenstates is lost be- 
tween successive resonances. For neutrinos that undergo 
a single H and a single L resonance, this is a valid as- 
sumption since the two resonance regions are well sepa- 
rated [3]. However, when the neutrinos encounter mul- 
tiple H resonances that are close to each other, coherence 
between the neutrino mass eigenstates is maintained, and 
one has to compute the amplitudes of neutrino flavor con- 
versions at the resonances, keeping track of the relative 
phases. This gives rise to the phase effects that we ex- 
plore in this paper. Such effects were pointed out in the 
context of solar neutrinos in and also observed in 
the numerical studies of supernova neutrinos in ^, ITsj . 
We provide an analytical approximation to study these 
effects, relating them to the the density profile of the 
medium. We also provide criteria to decide when deco- 
herence or finite energy resolution of detectors renders 
these effects unimportant. 

If the multiple resonances are semi-adiabatic, the phase 
effects may be strong and if the oscillations are indeed 
observed in the neutrino spectra, they can help us infer 
about the density profile of the shock wave, which in turn 
can provide us important clues about the shock wave 
propagation and the SN explosion. 

The observability of these effects is dependent on fac- 
tors like neutrino fluxes at the source, possible coher- 
ent development of neutrinos due to neutrino-neutrino 
interactions, the different density profiles encountered 
by neutrinos coming from different regions of the ncu- 
trinosphere, stochastic fluctuations of density, effective 
luminosity of the supernova, detector capabilities, et 
cetera. We comment on different model predictions of 
source fluxes and point out the model independent fea- 
tures of the phase effect. We investigate the number of 
events and detector capabilities required for the effect to 
be observable and find them to be demanding. 

Our emphasis in this paper is to first ascertain the 
origin of these phase effects, their dependence on var- 
ious parameters in the problem, and to check whether 
they are important at least in a simplified analysis where 
effects like coherent development of neutrinos [H, |T3], 
anisotropy f§j and turbulence [l^, of the medium are 
neglected. We find that even in a simplified setting, the 
observation of these effects is challenging. 

The paper is organised as follows. In Sec. |TT1 we give 
an analytical approximation for calculating the neutrino 
conversion probability when multiple H resonances are 
taken into account in a two-neutrino framework. The 
results of this section are general and can be applied to 
any situation where multiple resonances are involved. In 
Sec. mil we apply the results to the case of a SN shock 



wave, where the third neutrino and the L resonance are 
included. We show the neutrino conversion probabilities 
obtained by using a realistic shock wave profile and study 
the feasibility of detecting the phase effects. Sec. |IV] con- 
cludes. 



II. PHASE EFFECTS FROM MULTIPLE 
RESONANCES 

In this section, we calculate the survival probability of 
when they pass through multiple resonances, keeping 
track of the relative phases between the mass eigenstates. 
The calculations are performed in the 2-iy framework. 
The results are readily extended to the S-v framework in 
the case of a SN shock wave, as will be shown in Sec. IIIII 
Although all the arguments in this section are given for 
neutrinos, they are valid for antineutrinos just by chang- 
ing the sign of the matter potential A{x). However, the 
fiavor conversion analysis at the H resonance is appli- 
cable to neutrinos in the normal mass hierarchy and for 
antineutrinos in the inverted mass hierarchy. Therefore, 
while analyzing antineutrinos, the value of Am^ should 
be taken to be negative. 



A. Survival probability of for a small mixing 
angle 

The relevant mixing angle at the resonance H is 
9 = 9i3, on which we currently have a strong bound: 
sin^eiia < 0.05 T^. Therefore, we try to solve the prob- 
lem using a small angle approximation. We follow the 
notation and framework outlined in and work in the 
fiavor basis. 

Let be the relevant linear combination of v^, and v-r. 
When neutrinos pass through matter, their propagation 
is described by the Schrodinger equation 



dx 



H 



(2) 



Upto a matrix proportional to the unit matrix, the Hamil- 
tonian H is given by. 



1 / A{x) ~ Am'^ cos 29 Am^ gin 29 



AE 



Am^ sin 29 



'A{x) + Am^ cos 26* 



, (3) 



where A{x) = 2EV{x) = 2V2GFYepE/MN. These two 
coupled first order equations give rise to the second order 
equation 



(4) 



where 



1 Am^ 
(x) ^ —{A{x) - Am' COS29) , r^=^sm29 (5) 
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and prime (') denotes derivative with respect to x. In 
order to find the survival probabihty of v^, we solve for 
the Ve wavefunction with the initial conditions i^e(O) = 1, 
i'^(O) = 0. These conditions are equivalent to 



0(0). 



(6) 



The "logarithmic perturbation" approximation [20| 
solves the differential equation ^ for small mixing angles 
by choosing 

g = l-cos2e' (7) 
as the small expansion parameter. Denoting 

I., = e^(''), with S'{x) ^co{x)+gci{x) + 0{g^) (8) 
so that 

i^eix) = exp dxiCo(a;i) 

+ gj^ dxiciixi) + 0{g^)^ , (9) 



the solution becomes 

iQ{x) 



exp 



9- 



iAm X 



'2 {2E f Jo 



+0{g^) . (10) 
Here we have defined the "accumulated phase" 



1 /"^ 

Qix) = dxi (A(.Ti)- Ar 



(11) 



The survival probability Pee{x) = P^Ve — > t'e) at a; = X 
then becomes 



PeeiX) = exp 



(Am 



2\2 



2 {2Ey 



X 



dxie"^^""'' 



-0{g' 



(12) 



The integral in the above expression can be evaluated 
using the stationary phase approximation. The integral 
oscillates rapidly unless Q'{x) ~ 0. So the entire contri- 
bution to the integral can be taken to be from the saddle 
point Zs, which is the point where Q'{xs) = 0, i.e. 



A{xs) = Am^ 



(13) 



Note that this is also the resonance point in the small 
angle limit. 

For a monotonic density profile, there is only one sad- 
dle point Xg and the survival probability is 



exp 



7r(Am2)2 
^2E\A'{xs: 



(14) 



which agrees with the Landau-Zener jump probability 
pll [2^ in the limit of small mixing angle, and hence 
small (7, even when Pee ^ 1. 

For a non-monotonic density profile, neutrinos can ex- 
perience more than one resonance at the same density 
but at different positions. In that case Q'{x) = at 
more than one point. If the resonances are sufficiently 
far apart, the contributions from each of them may be 
added independently of each other. Their total contribu- 
tion to the integral in is 



X 



AttE 



1/2 



(15) 



where i runs over all the saddle points. Note that a = 
7r/4 if A'{xs) < and a = 37r/4 if A'{xs) > 0. The 
probability calculated using (fT2|) now has terms which 
depend on the differences between the integrated phases 



In general, 



Q{xj) - Q{xi) + < 



2E 



{A{x) - Am^) dx + aj - ai . (16) 



PeeiX) = exp 



-g j + 2 a-iaj cos 



where 



+0{g') 



ai = 



1/2 



2E\A'{x^ 



(17) 



(18) 



For example, when there are only two saddle points 
the survival probability is given by 



Pee = exp(— ga^) exp(— 502) exp(— 2(7aia2 cos $12) 



(19) 



The first two factors in are the individual Landau- 
Zener jump probabilities for the two level crossings. The 
last factor gives rise to oscillations in Pee as a function of 
energy. The oscillation pattern has its maxima at $12 = 
(2n-f l)7r and minima at $12 = 2mr where n is an integer. 
This expression is valid as long as gaf < 1, so that the 
0{g'^) terms remain small, and none of the resonances 
overlap. 

We illustrate the validity and limitations of the small 
angle approximation with a toy density profile 



p{x) 



a + hix"^ {x < 0) 
a + b2X^ (x > 0) 



(20) 



as shown in figure [TJ We take Ye = 0.5 and Am^ — 
0.002 eV^. Neutrinos are produced at x ^ —00 and we 
calculate Pee at a; ^ 00. We also show the positions of 



4 



o 
3 



0) 
Q 



6000 
5000 
4000 
3000 
2000 
1000 




^1 1 1 1 ; 1 

\ / P5 










Pio ■ 










P52 ■ 




-0.2 


-0.1 





0.1 


0.2 



Q. 
>. 

nJ 
o 

> 



1 

0.8 
0.6 
0.4 
0.2 






1 








1 







x(10° cm) 



10 20 30 40 
Energy (MeV) 



50 



60 



FIG. 1: Density profile in (|20p witii a = 500 g/cc, bi — FIG. 3: Survival probability Pee as a function of energy for 

10^ g/cc /(ID® cm)^ and 62 = 2.5x10^ g/cc /(ID* cm)^ . The 9 = 0.1 rad « 5.7°. The convention for the lines is the same 

horizontal lines on the graph are the resonance densities for as that used in figure (2] 
various energies (5, 10, 52 MeV) taking 8 = 0.02 rad 1.1°. 
Notice how the saddle points come closer for larger energies 
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FIG. 2: Survival probability Pee as a function of energy for 
= 0.02 rad ~ 1.1°. The green (light) curve is the numer- 
ically evaluated exact result. The blue (dark) curve is our 
solution with small angle approximation including the phase 
effects. The red (dotted) curve is the approximate solution if 
the phase effects are neglected. Notice that our approximate 
solution is valid only upto E = Ej^^^ax)- 



resonance densities for various energies, which are given 
by 



PB.[g/cc] « ± 



Am^[eV^] cos26> 
2 X 7.6 X 10-^YeE[MeY] 



(21) 



Figure [5] shows the survival probability Pee as a func- 
tion of energy, both the exact numerical result and the 
result of our analytical approximation for small angles. 
It can be seen that at such small angles {9 = 0.02 rad 
« 1.1°), the approximation works extremely well. 

Note that the amplitude of the oscillations is compa- 
rable to the deviation of the average survival probability 
from unity. That is, the oscillation effect is not a small 
effect. Indeed, the oscillation term is of the same order 
as the averaged effect, as can be seen from (fT9|) . Fig- 
ure [2] also shows the average value of Pee that one would 
have obtained if one naively combined the jump probabil- 
ities at the two resonances. Our analysis gives additional 



oscillations in the survival probability as a function of 
neutrino energy about this average value. This effect is 
what we call as the phase effect, and is clearly significant 
as can be seen from the figure. An important feature of 
the oscillations is that the "wavelengths," i.e. the dis- 
tances between the consecutive maxima or minima, are 
larger at larger E. 

The resonances start overlapping at £^ « ^'^(maa;) (fig- 
ure [T]), which is where our approximation starts breaking 
down, as can be seen in figure [51 For E > i?ij(maa;)j 
the neutrinos no longer encounter a strict resonance, and 
our approximation gives Pee = identically. However, 
the resonances have finite widths which may affect the 
conversion probabilities of neutrinos with E « Ej^j^ax)- 
The sharp jump observed in figure [2] at « 52 MeV is 
therefore not a real effect, but a limitation of our tech- 
nique. 

The small angle approximation starts failing for larger 
angles and lower energies, where gaf > 0(1). Figure [3] 
shows that the amplitude at low energies is not calcu- 
lated correctly for 6 — 0.1 rad « 5.7°. However, note 
that the positions of maxima and minima of Pee are still 
predicted to a good accuracy. We shall argue in the next 
subsection that these can be computed accurately for the 
whole allowed range of ^13, given the non-monotonic den- 
sity profile between the two resonances. 



B. The oscillations in Pee{E) 

Let us consider a density profile with a "dip" as in 
the toy model in the previous section. A neutrino with 
energy E encounters two resonances Ri and R2 at x = xi 
and X = X2 respectively, so that 



PR = p{xi) = p{x2) 



(22) 



We assume Ye to be a constant throughout the region of 
interest. We also assume that the propagation of neu- 
trino mass eigenstates is adiabatic everywhere except in 
the resonance regions (xi_,a::i+) and {x2-,X2+) around 
the resonance points xi and X2 respectively. In the limit 
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of a small mixing angle, the widths of the resonances are 
small: 



Ap « p tan 26* 



(23) 



Therefore, xi- « « xi+ and X2- ~ a;2 ~ X2+- We 
shall work in this approximation, and shall use the nota- 
tion Xi± only for the sake of clarity wherever needed. 

At a; <C xi, the density p{x) ^ pa, so that the heavier 
mass eigenstate is approximately equal to the flavor 
eigenstate Ve- Let us start with Ve as the initial state: 



i^eix < Xi) « l^H 



(24) 



The mass eigenstate i^h propagates adiabatically till it 
reaches the resonance region x xi: 



(25) 



While passing through the resonance, unless the reso- 
nance is completely adiabatic, the state becomes a linear 
combination of vh and j/l, the lighter mass eigenstate. 
Note that the phases of vh and vl can be defined to 
make their relative phase vanish at x = 2:1+. 



Ve{xi+) = cosxiJ/ff + sinxi^-L 



(26) 



where Pi = siv? xi is the "jump probability" at Ri if it 
were an isolated resonance. 

The two mass eigenstates uh and vl propagate to the 
other resonance i?2) gaining a relative phase in the pro- 
cess (the overall phase of the state is irrelevant): 



Vf,{x2-) = cosxit^ff + sinxi cxp 



2E 



dx j vl , 
(27) 

where Am^ is the mass squared difference between vh 
and z^L in matter: 

Ato^ = ((Am^ cos 29 - 2EV{x)f + (Am^ sin 261)) . 

The effect of the resonance R2 may be parametrised in 
general as 



yH[x2+)\ _ 
vl{x2+. 



cos X2 

' sin X2e 



-lip 



sm X2e^ 
cos X2 



vh{x2-) 
vl{x2-) 



(29) 

where P2 = sin^ X2 is the "jump probability" at R2 if it 
were an isolated resonance. 

From (jlOp . one can deduce that in the limit X2- ~ X2+, 
we have ip ~ Q{x2+ — X2-) ~ 0. The state h'e{x2+) can 
then be written as 



l^eix2+) 



cos X2 COS Xl 

'Sinx2 sinxicxp 

COS X2 sinxi exp 
smx2C0SXi VL . 



Arh^ 
2E 



-dx 



2E 



-dx 



(30) 



For x > X2+ , the mass eigenstates travel independently 
and over sufficiently large distances, decohere from one 
another. At x ^ X2, since p{x) 3> p_R, the heavier mass 
eigenstate vu again coincides with and we get the 
survival probability as 



Pee = COS^(xi - X2) 

— sin 2xi sin 2x2 sin^ 



Am2 



'iE 



dx . (31) 



If the phase information were lost, either due to deco- 
herence or due to finite energy resolution of the detectors 
(23! the survival probability would have been 

ne{no phase) = P1P2 + {I ~ Pl){l ~ P2) 

= cos^ Xl cos^ X2 + sin^ xi sin^ X2 , (32) 

which matches with (f3T|) when the sin^(/ ..) term is av- 
eraged out to 1/2. 

The sin^(/..) term in (|3ip gives rise to the oscilla- 
tions in Pee{E). If two consecutive maxima of P^e are 
at energies Ek and Ek+i such that Ek > Ek+i, then the 
condition 



Am''ix,Ek+i) , r^^^'-'^ Am\x,Ek) , 
dx — dx 



2tt 



2 Eh 



(33) 



is satisfied. The quantity {Ek—Ek+i) is the "wavelength" 
of the oscillations. 

Note that Afh?{x,E) is equal to \A{x) — Am?\ in the 
small angle limit. Moreover, this quantity is rather in- 
sensitive to 9 in the allowed range of 6'i3. Therefore, it 
is not a surprise that the predictions of the positions of 
maxima and minima in the small angle approximation 
(Sec. Ill Ap are accurate and robust in the whole range 
9 < 13°. 

Since 9 is small, the left hand side of (|33l) is approx- 
imately equal to the area of the region in the density 
profile plot enclosed by the densities pE^. and PEk+i- 



Ak,2-k 



M 



N 



V^GpYe 



(34) 



The distance between the two resonances in the region 
PSfc < P < PEk+i is then 



Tk ~ A/{pi 



PEk 



(35) 



This procedure may be repeated for various k values to 
estimate the separation between the resonances at the 
corresponding densities, and hence to constrain the form 
of the density profile. Although this seems straightfor- 
ward in our simplified analysis, the effect of density vari- 
ations due to convection, turbulence and anisotropics 
would greatly complicate such a reconstruction in prac- 
tice. 

Nonmonotonic density profiles are encountered by neu- 
trinos escaping from a core collapse supernova during the 
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shock wave propagation. If the phase effects are observ- 
able at neutrino detectors, the above procedure may help 
us reconstruct the shock wave partially. Of course the 
procedure is somewhat crude, and the information ob- 
tained would only be on the density profile along the 
line of sight. We have also assumed that neutrinos com- 
ing from different parts of the neutrinosphere encounter 
nearly the same density profiles. However, since neu- 
trinos are the only particles that can even in principle 
carry information from so deep inside the exploding star, 
it is important to check whether the detection of these 
phase effects is feasible even with these simplifying as- 
sumptions. We shall do this in the next section. 

III. OSCILLATIONS DURING THE SN SHOCK 
WAVE PROPAGATION 

In this section, we apply the results in the last section 
to the neutrinos travelling through a supernova shock 
wave. Though we have to consider three-neutrino mix- 
ing in this case, the separation of H and L resonances 
[24| means that we can calculate the transition probabil- 
ities at these resonances separately. Each of these reso- 
nances can then be treated as an effective two-neutrino 
level crossing. The L resonance that takes place in neutri- 
nos is always adiabatic ^] . If the net survival probability 
after passing through all the H resonances is denoted by 
Ph, the survival probability of i/g after passing through 
all the H and L resonances is 

p = P//sin2 6'o (NH), p = sin^ 6*0 (IH) (36) 

where NH and IH stand for normal and inverted mass 
hierarchy respectively. Here 9q is the solar mixing an- 
gle. Similarly, the survival probability of Ve after passing 
through all the H and L resonances is 

P=Phcos^Oq (IH), P^cos^Oq (NH) . (37) 

Clearly, since the phase effects appear through Ph, they 
will be visible only in Ve for normal hierarchy and only 
in De for inverted hierarchy. 

In present and planned water Cherenkov [25| and scin- 
tillation [2^ detectors, the main neutrino detection chan- 
nel is the inverse beta decay reaction Vf,p — > ne"*" that 
allows the reconstruction of v^. energies. Therefore we 
consider only the Vf. spectrum in our analysis. However 
an analogous analysis can be easily performed in the neu- 
trino channel for a detector able to measure the spec- 
trum, for example using liquid argon [23|. 

In order to illustrate the phase effects on p or p, we 
consider a typical snapshot of the density profile of a SN 
during a shock wave as shown in figured) The for- 
ward shock F and the reverse shock R are sharp density 
discontinuities, the density change of a factor of two or 
more taking place over a distance of much less than a 
km. The density variation in the "contact discontinu- 
ity" C, which is the transition region between the shock- 
accelerated and neutrino-heated SN ejecta, takes place 
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FIG. 4: Snapshot of a shock wave density profile at t = 5 
sec. The resonance densities for E = ^ and 80 MeV (with 
Ye = 0.5) are also shown. We take Am^ = 0.002 eV^ and 
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more slowly, over a distance of more than 100 km [8| . The 
mass accretion region A behind the forward shock wave, 
and the low density bubble B have gradually changing 
densities. The region T is the tail of the shock wave. 

The neutrinos, while passing through these regions, 
may undergo multiple level crossings. The extent of fla- 
vor conversion in each region will depend on the value of 
^13 and the steepness of the density profile in that region. 
It is found that for ~ 0.01 rad or higher, the density 
variations in the mass accretion region A, the low density 
bubble B and the contact discontinuity C are too grad- 
ual for any non-adiabaticity. We therefore concentrate 
on the forward shock F, the reverse shock R and the tail 
T. 

The coherence between mass eigenstates and the oscil- 
lations in the survival probability may be lost due to 
two major sources: the separation of mass eigenstate 
wavepackets and the finite energy resolution of the detec- 
tor. The coherence length Lcoh, defined as the distance 
over which the wavepackets separate, is given by [1^ 

Lcoh - 4V2CTi;VAm^ , (38) 

where a is the width of the wavepacket at source. Taking 
a ^ 10~^ cm near the neutrinosphere p9l | in the relevant 
energy range of 5-80 MeV, the coherence length for SN 
neutrinos is Lcoh ~ lO'^-lO^ km. Resonances separated 
by distances larger than Lcoh may be taken to be inco- 
herent. Since the distances involved are 0(10'' km) (See 
figure S]), a definite conclusion about decoherence due to 
wavepacket separation cannot be reached with this sim- 
ple estimate. However, for observability the oscillations 
in survival probability must also occur over energy inter- 
vals much larger than SEdet, the uncertainity in energy 
measurement at the detector. This turns out to be the 
dominant factor in smearing out the oscillations. We can 
estimate from the difference in energies at which 
successive maxima of the survival probability occur to 
be 

Xe, = Ek+i ~Ek^ 4TTEl/Am\k , (39) 
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where is the distance between the two resonances en- 
countered by a neutrino of energy E]^. For the energy 
range of 5-80 MeV, this gives A^^ of 1-10 MeV for the 
T-R resonance pair with ^ lO'^-lO^ km. The oscil- 
lations are faster for the T-F and R-F pairs that occur 
about 10^-10^ km apart. Typically the energy resolution 
is 1-10 MeV for a water Cherenkov detector and 0.1-1 
MeV for a scintillation detector over this energy range. 
Moreover, the charged lepton energy is not of the same as 
the energy of the neutrino, which introduces additional 
smearing. Indeed as we shall see, the fast oscillations 
due to the T-F or R-F resonance pairs will be smeared 
and therefore clearly unobservable even at a scintillation 
detector. Note that the above arguments are only qual- 
itative; we have neglected the matter effects while esti- 
mating hcoh and ■ 

Figure [5] shows the value of p as a function of energy for 
= 0.02 rad k, 1.1°. The rapid oscillations correspond to 
the relative phase ^rf{~ ^tf) that is accumulated by 
the mass eigenstates between resonance regions R and F 
(T and F). Such high frequency oscillations are virtually 
impossible to observe, given the practical limits on the 
energy resolutions of neutrino detectors. 

We smear the high frequency oscillations by taking a 
"running average" over the energy range corresponding 
to the typical energy resolution of a scintillation detector. 
The low frequency oscillations that survive are found to 
correspond to the relative phase '^tr accumulated by the 
mass eigenstates between resonances in regions T and R. 
Since these two resonances are closer compared to the 
resonance pairs R-F or T-F, the frequency of oscillations 
is smaller. The same oscillation pattern is observed if 
the survival probability is computed by assuming that 
the resonance in region F is completely non-adiabatic, 
which confirms that the pattern is indeed due to the level 
crossings in regions T and R. 

Note that p goes to cos^ 9q at its maximum where Pr 
goes to unity, whereas at p « at high energies and low 
energies where Ph goes to zero. The oscillations in the 
low energy region [E < 20 MeV) are too rapid to be 
observable. The fluctuations observed in the running av- 
erage at low energies are not robust: they depend partly 
on the details of the density profile and are partly nu- 
merical artifacts. 

In figure [SI we plot the value of p for the same pa- 
rameter values, but calculated using ([TT]) at the T, R, F 
resonances. The smeared probability curve is calculated 
by dropping the oscillatory terms due to the resonance 
at F and combining the survival probability at the T-R 
pair with that at F using ([5^ . We see that the ana- 
lytical expression agrees quite well with the numerical 
result of figure El In particular, the oscillation frequency 
matches quite well. Moreover, the slope of A{x) at T is 
about 0.1 times that at R, therefore (US]) and ([TO]) predict 
the amplitude of oscillations to be ^y A' (T) / A' {R) « 0.3 
times the mean value of p. This estimate is also in good 
agreement with our numerical results in the figure. 

The analytical approximation in (I17p breaks down 
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FIG. 5: Survival probability p of i>e for inverted hierarchy 
with the density profile in figure U The smeared probability 
is obtained by taking a running average over the energy range 
corresponding to the typical energy resolution of a scintilla- 
tion detector i.e., A£;sc(MeV) « 0.2y/E/W MeV. 
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FIG. 6: Same as figure O calculated analytically using (|17|) 
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dropping the fast oscillatory terms. Note that analytical ap- 
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FIG. 7: Comparison of smeared probability calculated nu- 
merically and using the analytical approximation, for a larger 
mixing angle 613 = 0.05 rad « 2.9°. Note that analytical ap- 
proximation reproduces the numerical result reasonably well 
for E >25 MeV. 



when gaf > 0(1). This happens below a certain value of 
energy, that is higher for larger mixing angles. Compar- 
ing figure [5] and figure [S] we see that for 9 — 0.01 rad the 
approximation is resonably accurate for neutrino energies 
above 20 MeV. Figure [7] shows the smeared p calculated 
numerically and analytically for 9 = 0.05 rad. It is clear 
that the approximation is valid above 25 MeV. The ap- 
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FIG. 8: The upper figure (a) shows the osciUations in the 
energy spectra of the Ue, if the Ve energy were accurately 
measurable, for the Garching and Livermore models. The 
bottom figure (b) shows the observed positron energy spectra 
in a scintillation detector. All the spectra are normalized 
so that the total number of events is 10^. The figure (b) 
also shows the statistical fluctuations for 10^ events with an 
energy binning of 0.5 MeV. 



proximation fails at low energies where p calculated an- 
alytically does not tend to zero. For higher energies, the 
oscillations due to the phase effect are predicted quite 
accurately. 

In order to explore the observability of the phase ef- 
fects, we use the parametrization for the primary fluxes 
given by [sot : 



exp 



-(" + i)7r 



(40) 



where A^(a) = (1 + a)i+"/r(l + a). For illustration, 
we choose two models of neutrino fluxes, the Garching 
model [3l| that uses the parameters 

ape ~ Q!px = 3 , Eq{v^) = 15 MeV , 
Eo{V^) = 18 MeV , M^e)/M^x) = 0.8 , 

and the Livermore model 13211 that uses 



ap, = ap^ = 3 , EQ{ve) = 15 MeV , 
E^{i?^) = 24 MeV , $o(^>e)/$o(^^.) = 1-6 . 

We plot in figure [H^a) the quantity 



(41) 



= a{E,^)[Fl+p{Fl-Fl)] (42) 



as a function of the neutrino energy E^^ , where we have 
normalized the spectrum such that the total number of 
events is 10^. Here we use the differential cross section 
as computed in [33]. 

The phase effect is most prominent around 4-5 seconds 
postbounce when the tail T and reverse shock R pass 
through the H resonance and the T-R resonance pair 
is about 10'^ km apart. We expect from ([55)1 that the 
extrema in the survival probability will occur at energies 
separated by about 1-10 MeV. This is clearly visible in 
figure [5fa). Moreover, the positions of the extrema are 
independent of the primary fluxes. 

However, the spectrum shown in figure [HKa) is not di- 
rectly observable: one can only observe the energy spec- 
trum of positrons produced by the inverse beta reaction. 
Assumin g qu asielastic scattering, the positron energy is 
given by [33| 



_ {Eu + + ecos^^{E, - 5Y - mln ^^^^ 



and 



Pi 



= Ve! 



(44) 



where is the angle of scattering, Mp the proton mass, 
e = E^JMp, and k = {I + - {ecos^f. For E k 4Q 
MeV, we have e w 1/25, so that the positron energy 
is spread over a range of ~ 4 MeV depending on the 
scattering angle. Given that the successive maxima of 
the oscillation pattern are separated by only about 2- 
8 MeV in this energy range, the oscillation pattern is 
significantly smeared out. 

A further smearing of the oscillation pattern is caused 
by the finite energy resolution of the detector. The en- 
ergy resolution of a water Cherenkov detector is typically 
A£;cH(MeV) « l.&yfE/lQ MeV and washes off the os- 
cillations completely. For a scintillation detec tor, the res- 
olution is much better, A£'sc(MeV) w 0.2^£:/10 MeV. 
We show in figure [5{b) the spectrum of the observed 
positron energy E'^ after taking and ([33]) into ac- 
count and using the energy resolution of a scintillation 
detector. It is observed that one or two extrema at high 
energies {E k, 40-60 MeV) may still survive for the Liv- 
ermore model where the spectrum extends to higher en- 
ergies, but their clean identification would require ~ 10^ 
events at a scintillation detector in a single time bin. 
Figure ^ih) also shows the positron energy spectrum for 
10^ events, binned in 0.5 MeV energy intervals, which 
is approximately the energy resolution of a scintillation 
detector near E = MeV. 

The total number of events expected even at a future 
50 kt scintillation detector is ^ 10^ for a SN at a distance 
of 10 kpc f26j. Even with 1 sec time bins, the number 
of events in each bin will be ~ 10''. This is a number 
too small for our purposes. Thus, the identification of 
the phase effects seems very unlikely, unless the SN is as 
close as a kpc. 

The density profiles in the shock wave are quite un- 
certain, and one may expect that some possible profiles 
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give rise to oscillations with larger amplitudes and wave- 
lengths, which would be easier to observe. The amplitude 
of oscillations is proportional to the ratio \ A'{T)/A'{R)\ 
[See (dH]) and HH)]. Therefore, larger amplitudes need 
a larger ratio \ A'{T)/A'{R)\ whereas larger wavelengths 
need a sharp tail, i.e. a larger |A'(r)|. At the same time, 
the adiabaticity parameter at the two resonances has to 
be 0.1 < 7 < 2 for E 40 MeV. However, we find that 
even with such finely tuned density profiles, the improve- 
ment in observability is not significant. 

IV. CONCLUDING REMARKS 

When the neutrinos escaping from the core of a core 
collapse SN pass through the shock wave, they may 
encounter multiple "i?" resonances corresponding to 
Am^^jjj and ^13, when the shock wave is in the regions 
with densities around 500-5000 g/cc. We have shown 
that this necessarily gives rise to oscillations in the neu- 
trino survival probabilities, which we have calculated as 
a function of energy. We present an analytical approxi- 
mation for small mixing angles and show that the oscilla- 
tions are a significant effect: they can be of the same or- 
der as the non-oscillating terms. The typical values of ^13 
that gives rise to the oscillation features are in the "tran- 
sition region" of the neutrino mixing parameter space Q 
that is significant in range (10~^ < sin^ 813 < 10"'^) but 
is usually neglected in the SN analysis for the sake of 
simplicity. 

The local maxima and minima in the survival proba- 
bilities of Ve or i>e are determined by the relative phase 
accumulated by the neutrino mass eigenstates between 
the multiple resonances. The positions of these extrema 
depend on the density profile and are independent of the 
primary neutrino spectra. If these extrema were identi- 
fied they would reveal information on the propagation of 
the shock wave: its location as well as the density varia- 
tion present around it. Since neutrinos are the only par- 
ticles that can carry information about the shock wave 
while it is still deep inside the exploding star, it is im- 
portant to explore the observability of these phase effects. 
Moreover, the mere identification of these effects in the 

spectrum would establish the inverted hierarchy and 
nonzero ^13, which are two of the most important quan- 
tities in neutrino phenomenology, and even faint chances 
of their determination should be explored. 

It is interesting that for typical shock wave profiles, 
oscillations with "wavelengths" of 2-8 MeV are indeed 
present in the neutrino survival probabilities in the en- 
ergy range E « 30-60 MeV. These wavelengths are tan- 
talizingly close to the resolving power of the neutrino 
detectors. However, the inability of the detectors in re- 
constructing the energy of the incoming neutrino tends to 
wash out the oscillation pattern. For typical shock wave 
density profiles, the energy resolution of water Cherenkov 
detectors is insufficient to detect the oscillations. Even 



a scintillation detector with a superior energy resolution 
will need ~ 10^ events in the relevant time bin of 4-5 
sec postbounce for identifying one or two extrema in the 
most optimistic scenario. 

Therefore, we expect that neglecting the phase effects 
is a valid practical approximation except under extreme 
cases. Note that it is still possible to observe other robust 
signatures of shock wave propagation like the dips in the 
time spectrum of number of events Q or the double dip 
feature in (E) Q. Some of these features survive even 
in the presence of somewhat extreme stochastic density 
fluctuations, as discussed in [l^]. The oscillatory effects 
might contribute additional scatter to these signals but 
will not spoil the signatures. 

We must qualify the above conclusions by stating the 
effects ignored in our analysis to obtain the present re- 
sults. We have assumed a smooth spherically symmetric 
density profile and ignored anisotropics that are likely 
to be present [^]. This is a justified assumption only as 
long as the deviations from this assumed profile occur 
only on transverse distances larger than the size of the 
neutrinosphere so that the neutrinos in our line of sight 
do not experience the anisotropy. We have not included 
the recently discussed collective effects of coherent fla- 
vor development [T7| , which may be important in the 
inverted hierarchy. However, the extent and nature of 
its impact on flavor conversion has not been worked out 
in detail, which makes it difficult to be included in the 
present analysis. Similarly, the effects of a realistic spec- 
trum of stochastic fluctuations in the medium density or 
turbulent convections behind the shock wave are yet to 
be calculated for SN neutrinos [l^, [ll|. We think that 
with a better understanding of the collective effects and 
density fluctuations one could include their effect on the 
realistic observability of the phase effects. 

The phase effects pointed out here result from the in- 
terference between two or more MSW resonances. This 
phenomenon is not restricted to the SN alone, but may 
occur whenever neutrinos pass through nonmonotonic 
matter densities and the resonances are semiadiabatic. 
The technique developed in this paper for treating co- 
herent multiple resonances is applicable to such cases. 
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